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Abstract We present a theory for Coulomb drag between 
two mesoscopic systems which expresses the drag in terms of 
scattering matrices and wave functions. The formalism can 
be applied to both ballistic and disordered systems and the 
consequences can be studied either by numerical simulations 
or analytic means such as perturbation theory or random 
matrix theory. The physics of Coulomb drag in the meso- 
scopic regime is very different from Coulomb drag between 
extended electron systems. In the mesoscopic regime we in 
general find fluctuations of the drag comparable to the mean 
value. Examples are the vanishing average drag for chaotic 
2D-systems and the dominating fluctuations of drag between 
quasi-ballistic wires with almost ideal transmission. 

1 Introduction 

Current flow in a conductor can through a Coulomb me- 
diated drag-force accelerate charge-carriers in a nearby 
conductor, thus inducing a drag-current. The effect is 
active whenever the distance between the two conduc- 
tors is of the same order as the distance between the 
charge-carriers - otherwise it is suppressed by screen- 
ing. In the past years Coulomb drag in extended 2D- 
systems has been studied extensively and very re- 
cently the study of fluctuations of the Coulomb drag 
has been initiated As in many other mesoscopic 

phenomena the fluctuations will be pronounced for tem- 
peratures smaller than the Thouless energy. We study 
drag in disordered or chaotic finite-size systems with di- 
mension much smaller than the phase-coherence length 
iij) (see Fig. 1) and find interesting phenomena like sign- 
reversal, vanishing mean value, and large fluctuations of 
the drag response. 

2 Formalism 

Starting from the Kubo formula (as in Refs. [^,D, but, 
here we consider systems with broken translation invari- 
ance as illustrated in Fig. 1) we calculate the drag con- 
ductance G21 to second order in the interaction U12 be- 
tween mesoscopic subsystems, taking the isolated sys- 




tems to be otherwise non-interacting. In the dc limit 
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Here, 



A,{uj, r, r') = -2nr^h ^ 6'^(r, r', - hto) 

X [npi^ep - huj) - nF(£/3)] + (r ^ r'; a; -lu), (2) 
is the three point correlation function (Ipp) and 
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where i labels the subsystem and ^q, = The matrix 

elements are /^^ = {a\P\j) and Pafj{r) — (a|r)(r|/3), 
where |q;)'s are the eigenstates of the uncoupled subsys- 
tem with energies . With scattering states as the basis 
(r"^ — SW^S)^g is expressed in terms of 
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Fig. 1 Schematic geometry of a mesoscopic Coulomb drag 
experiment. 



the 2N X 2N scattering matrix S and the third Pauli 
matrix . 

Many interesting features follow immediately from 
eq. (|^) ||] but here we focus on two of them: 

— Phase-coherencen and lack of translation invariance 
leads to an arbitrary sign of the outcome of eq. (Q). 
This indicates that fluctuations might be pronounced 
in disordered or chaotic systems. 

— For T <^Tf the factor sinh^^ cuts off the frequency 
integration and Z\ (x w to lowest order. Performing 
the frequency integration gives G21 cx T^, in accor- 
dance with the usual Fermi liquid result for electron- 
electron scattering. 

In the following we study eq. (P in the T^-regime for 
two cases which demonstrate fluctuations of the order of 
or larger than the mean value. 

3 Quasi-ballistic one-dimensional wires 

We consider two weakly disordered ID-wires of sepa- 
ration d, length L <ti and mean free path £ ^ L. 
We study this system both analytically using perturba- 
tion theory and numerically by mapping the problem 
on to a lattice (see inset of Fig. ||). We flnd that even 
a very small amount of disorder {L <^ £) can lead to 
large fluctuations for the drag response and even re- 
verse the sign. The origin is inter-wire interaction in- 
duced forward scattering which gives rise to a drag re- 
sponse provided it is combined with disorder induced 
backscattering. In contrast, in the case of clean wires 
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Fig. 2 Panel (a) shows fluctuations as a function of kpi for kpL — 100 and kpd = 1. The fuU hues are Eqs. (^) and 
and the squares (triangles) are the numerically solution on a lattice (inset) for mutually uncorrelated (correlated) disorder. 
Panels (b) and (c) are the numerically obtained distributions (histograms based on ~ 10* random disorder configurations) for 
kpi ~ 3.8 X 10^ in the case of mutually uncorrelated and correlated disorder, respectively. Panel (d) shows the topologically 
different diagrams contributing to the fluctuations to lowest order in l/kpl- 



the backscattering is induced solely by the interwire in- 
teraction, and therefore the disordered case is larger by 
a factor of order (7^)[/l2(0)/C/l2(2/cF), with f7i2(g) = 

So Jo '^■^1 e^''^'^^~^^'^Ui2{xi,X2) being the Fourier trans- 
formed interaction. The reflection coefficient is inversely 
proportional to the mean free path (7^) ~ L/£ for L < £. 

Assuming that the disorder potentials of the two wires 
are mutually uncorrelated (uc) |^,||,D this can be shown 
explicitly by lowest order perturbation theory in l/kp£, 
corresponding to the diagram (3) shown in panel (d) of 
Fig. ^ (diagram (1) gives a vanishing contribution and 
diagrams (2) and (4) are not relevant to uncorrelated 
case). For kpL ^ 1 the result is 

2\l/2 r„/„ w„ w,2 toL_\rT2 /nM 1/2 



with Ui2{0) replaced by 
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The denominator is the result 6*21(00) cx Uf2{2kp) for 
ballistic wires. For the realistic case Uuiikp) ^ UuiO) 
it then follows that the fluctuations of the drag can ex- 
ceed the average value (G2i(£))^^ ~ G2i(oo). This con- 
trasts the behavior of the diagonal conductance where 

1/2 

the fluctuations (^[SGa]'^') are vanishing compared to 



the mean value (Gu) = {2e'^/h){l - (7^,)) - 2e'^/h in 
the limit of weak disorder. 

The opposite limit with mutually fully correlated (c) 
disorder potentials [^[sj is at first sight more compli- 
cated because the calculation of fluctuations involves 
also "crossed" diagrams like diagrams (2) and (4). How- 
ever, for correlated disorder both of the diagrams (3) 
and (4) contribute equally whereas for uncorrelated dis- 
order only diagram (3) is relevant. The same is the case 
for diagrams (1) and (2) and more generally, for each 
topologically different diagram contributing in the case 
of uncorrelated disorder there are two similar diagrams 
contributing equally in case of correlated disorder (to 
lowest order in l/kpt). Since symbolically (2) = (1) and 
(4) = (3) we get 



1/2 



^ (l) + (2)-K3) + (4) 
(1) + (3) 



V2. (6) 



{[^G2l{l)f) 

{[5G2i{lWfi^ 

Such an enhancement by correlated disorder was recently 
predicted for the mean drag in 2D systems 0; here we 
find that also the fluctuations are enhanced. 

Panel (a) of Fig. shows a parameter-free compari- 
son of Eqs. (^) and (^ to numerical results for kpL = 
100 and kpd = 1 in the case of bare Coulomb interac- 
tion. Both the l/fc^^£-dependence, the magnitude, and 
the \/2- enhancement are fully confirmed by the numeri- 
cal simulations and similar agreement has been found for 
other values of d and i . Panels (b) and (c) show typ- 
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ical numerically obtained distributions of G21 (^) for un- 
correlated and correlated disorder, respectively {kp£ ^ 
3.8 X 10^). In panel (b) the negative values correspond 
to sign-reversal of the drag. In panel (c) we see an en- 
hancement of both the mean value and the fluctuations 
compared to panel (b). Predictions for the full distri- 
butions do not yet exist and especially the absence of 
sign-reversal in panel (c) is interesting. For only partially 
correlated disorder we numerically find enhancements of 
the fluctuations in the range [1; -\/2] with the two limits 
corresponding to uncorrelated and fully correlated dis- 
order, respectively 

4 Chaotic mesoscopic systems 

We consider an ensemble of mesoscopic chaotic systems, 
such as suggested in Fig. |l|. We assume that the region 
where the subsystems couple by Coulomb interactions 
have mutually uncorrelated disorder so that (suppressing 
the integration variables) 



(G21) cx J Ui2Ui2{Ai){A2), (7) 
(G^i) cc J Ui2Ui2Ui2Ui2{AiAi){A2A2). (8) 

Starting from eq. (|^) in the low temperature limit the 
task is to calculate (Ai) and (Z\iZ\i). We do this using 
random matrix theory |9(] where the eigenvalues and the 
wave functions are uncorrelated. 

To lowest order in l/kp£ the average of two wave 
functions is (see e.g. Il0| ) 

{^;{x)My))^S^A{A{x-y)), (9) 

where Se is the level spacing and the spectral function is 
given by 

(A(r)) ~ (m/2n2)exp(-r/2^)Jo(ftFr). (10) 
Next, consider the average 

{Iai3(t)*ix)<j)s{y)) CX {d^, ~ d^^) (11) 
x{(t>l{xi)(f)p{x2)<f)*^{x)(t>s{y)) 
where to lowest order in l/kpl 

{<l>l{xi)M^2W^{^)^s{y)) ^ {rAxi)M^2)) (12) 

x{r,{x)^s{y)) + {rM^s{y)){Mx2W,{^)). 

Due to current conservation the points xi = X2 can be 
anywhere. Taking them to be outside the chaotic re- 
gion, the decay of the spectral function makes the sec- 
ond term in eq. ( p^ ) vanish if x or y are inside. The 
first term amounts to performing the average over Ia/3 
and (j)*{x)(l)s{y) separately. Similar arguments hold for 
higher order averages so that to lowest order in l/kp£ 
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and since / and p are Hermitian 



(13) 



{A,{uj,r,r')A,{Cj,s,s')) ^ 

^ ' a/37 5/37 

x{(4.4^)(pL/3(r)p^,(r-')p^^(^)p)^^(^')'5(^'')) 

-(^;,4«)(^'«/3(^)^^7(^')p^5(^)p;^(«')'5(e^)) 
+(^;.^;<i)(p^.(r-)p;^(r')p^<,(^)p;3(«')'5(c^))},(i4) 

where ^(e^) ^ <5(f^)<5(^J)<5(ef )<5(^|). 

From the statistical properties of the 5-matrix ^ we 
find {la-y} oc r^^ and since the second average in (A) 
is symmetric with respect to interchange of a and 7 we 
get {A) = and thus of course (G2i)=0. However, the 
fluctuations are nonzero. Since 

(lajla^i) = const. X T^^T?^ 

+ {4^f{{S^r'SUiS^r'Sh^), (15) 

and the average in the limit of a large N becomes 
(2iV)-i,5„^(557 we get |ll| 

{A,{uj,r,r')A,{iJj,s,s')) 1 



hujUoj {2Trffi^ ~ 4(27r)2?i2(2i)2(2iV)(27r)4 

x{{Mr,s')Mr\r)Ms,r')A{s' ,s)) 

-{Mr,s)Mr',r)Ms,s')Ms',r')) 
-{A{r, r')Mr', s')Ms, r)Ms', s)) 

+{Mr, r')Mr', s)A,{s, s')Ms' , r))}. (16) 

Here, we have introduced the spectral function A{r, r') = 
2TTY.o,(t>l{r)(t>a{r')5{£,^). To lowest order in l/kpl we 
replace each spectral function by its average and use 
that {A{r,r')) = {A{r' ,r)) so that 



{Ai{u,ry)A,{Co,s,s')) ^ T,{ry,s,s') 
■K^hujhCo ~ 2(2iV)(27r)4 

T.{ry,s,s') = {A,{ry)){A,{s,s')), 



(17) 
(18) 



(A,(r, s)){A,y, s')) ~ (A,(r, s')){A,{r' , s)) 



Since {{5G2if) = (G^i) - {G2if = {Gl^) we perform 
the (jj-integration in eq. (^ and get 

/ ,2x1/2 ^, {kTf 

Ui2{ri,r2)Ui2{r[,r'2)Ui2{si,S2)Ui2{s[,s'2) 

1/2 



xTi{ri,r[, Si,s[)T2{s2, s'^, r2,r'2] 



(19) 



To obtain an estimate we note that eq. (M) is very 
peaked for kpr < 1 and on long length scales^yl(r)) « 
(tt/Asp kpt)5{r) for kpt ^ 1. Using that approximation 
for all spectral functions is too crude, but 
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{A{ry)y 



X [5{r ~ s)5{r' - s') - 5[r - s')5{r' - s)] , (20) 
still gives a finite answer 



2 7 
1/2 TT^ 



kT 



h ?,2^N [kpiY \eF 



(21) 



{A{r^,r\))\A{s^,s',)y 



X [C/i2(ri, S2)[/i2(rl, s^,) - [/i2(rl, S2)C/i2(ri, s^)] ' 

Due to the peaked behavior of the spectral functions 
the mixed terms gives a vanishing contribution and in- 
troducing u = ri — r'l, V = S2 — and y = r[ — s'2 



1/2 



((^G2r)^)^/^-^' 



h 32«N{kFi)^ \eF 



kT 



2A 



X / {A{u,o)y{A{v,o)yuUu + y,v)uUy,o) 



(22) 
1/2 



where A is the interaction area. Assuming that the screen- 
ing length rg ^ i, the range of (A) is longer than the 
range of U12 and we further approximate 



((^G2l)^) 



2\l/2 



kT 



h i2^N {kpif \eF 



I2A 



di;(A(-u,0)) / AuUl2{u,^) \ dyUUy,0) 



(23) 

1/2 



The integral over the spectral function gives 

dr(^(r,0))' = 2. rdrr{A{r)y^^^-^^4^, 
Jo ^ £p 

where (^(kp^) — j'^ dxxe^'^^^^''^ Jq{x) can be approxi- 
mated by Cikpf) « (3/27r2)ln2fcF^ for kpi > 1. Simi- 
larly, for a screened interaction f 12 (r) = Ui2{r)exp{—r/rs), 
£^12(0 = e^/(47reoerr) we get 

J drUf2{r,0) = Stt ^ dr rUl2{V + d^) 
= 27rUUd)d^roi2d/r,), 

where ra{x) = dtt'^~^e~'^ is the incomplete Gamma 
function. Collecting things we now obtain the estimate 



,2\i/2 fkT Ui2{d) 



h \6p £p 

ro(ff)d2 kpVA \n2kpl 



(24) 



where 7 = 2'^^ ^/^^ ~ 0.7 x 10"^ |T|. Two sepa- 
rately contacted and electrostatically coupled quantum 
dots have recently been realized in a double-quantum- 
weh GaAs/AlGaAs with d 40nm and \fA 0.8 ^m 



13 . Transport measurements indicated a strong elec- 
trostatic interaction between the dots, though, measure- 
ments of Coulomb drag where not performed. For typical 
numbers we estimate the drag fluctuations to be of the 
order of 0.1 Ohm. 

Interestingly, the above results for the mean value 
and fluctuations are not changed by breaking of time 
reversal symmetry, in contrast to the UCF case, where 
the results with or without an applied i3-field differ by 
a factor of 2 |^ . 

5 Conclusion 

We have developed a formalism for study of Coulomb 
drag in the mesoscopic regime which is a promising new 
direction for the study of mesoscopic transport proper- 
ties, since it gives an opportunity to directly study inter- 
action and correlation effects in mesoscopic structures. 
Disordered mesoscopic systems exhibit interesting and 
unusual physics and we have demonstrated that this 
is also the case in disordered Coulomb drag systems. 

For quasi-ballistic ID-wires [L <C t) our results illus- 
trate how the statistics of the drag conductance depend 
strongly on disorder and we find that even weak disorder 
can give rise to fluctuations of the same order of magni- 
tude as the drag conductance for the ballistic case. This 
implies that the direction of drag depends on the disor- 
der configuration and that for a given system the sign 
of the drag current will be arbitrary. Depending on its 
nature the presence of disorder may even enhance the 
average of the drag response and the fluctuations com- 
pared to absence of disorder. 

For chaotic 2D systems like coupled quantum dots we 
find a vanishing mean drag but with finite fluctuations. 
These results are based on standard random matrix the- 
ory and the prediction of zero mean drag can thus be 
taken as a test of the degree of ergodicity of the system 
under investigation. 
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